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The spontaneous emission rate of a radiating atom reaches its time-independent equi- 
librium value after an initial transient regime. In this paper we consider the associated 
relaxation effects of the spontaneous decay rate of atoms in dispersive and absorbing 
dielectric media for atomic transition frequencies near material resonances. A quantum 
mechanical description of such media is furnished by a damped-polariton model, in which 
absorption is taken into account through coupling to a bath. We show how all field and 
matter operators in this theory can be expressed in terms of the bath operators at an 
initial time. The consistency of these solutions for the field and matter operators are 
found to depend on the validity of certain velocity sum rules. The transient effects in the 
spontaneous decay rate are studied with the help of several specific models for the dielec- 
tric constant, which are shown to follow from the general theory by adopting particular 
forms of the bath coupling constant. 
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I. INTRODUCTION 

The rate and the spectral and spatial characteristics of the spontaneous decay of an atom depend on 
the properties of the atom and of the radiation field, and on the interaction between them. The radiation 
field changes by the presence of other matter [Q . One can try and manipulate the emission properties 
once the influence of this medium is understood. 

In quantum optics of linear dielectrics, one tries to describe the material medium in an effective way 
with the help of the classical dielectric function e(r,oj), which in general is a complex function of both 
position and frequency and in this full generality describes the propagation and loss of light at each point 
in the dielectric. Sometimes it is possible to neglect the spatial variations (including local field effects), 
dispersion and losses altogether. The spontaneous emission rate of an atom in such a simple dielectric is 
the refractive index n of the medium times the rate Tq in vacuum |^| — Q|. 

The situation becomes more complicated when material dispersion has to be taken into account [|5j - 
PH| . Since the Kramers-Kronig relations tell that dispersion and loss always come together (be it not 
always at the same frequencies) , one should like to include losses as well in order to describe all frequencies 
in one theory. The damped polariton model jl2| - provides us with such a microscopic theory. From 
that theory it was shown that the radiative spontaneous emission rate equals Tq times the real part of 
the refractive index at the transition frequency Jl6| ]. 

The quantum mechanical treatment of dissipative systems is more complicated than the classical one, 
because of the extra requirement that equal-time commutation relations do not change over time [jr^.|l8| . 
Based on the damped-polariton model and on the fluctuation-dissipation theorem, phenomenological 
quantization theories were constructed that meet these requirements. In these theories, the dielectric 
function is an input function and the Maxwell field operators satisfy quantum Langevin equations with 
both loss and quantum noise terms p9pC| ] . With the use of a Green- function approach, the phenomenolog- 
ical quantization theories have been generalized to inhomogeneous dielectrics, first for multilayer systems 
and later for general e(r,u>) [|TJ - |£2). Field commutation relations turn out only to depend on the 
analytical properties of the Green function. However, the calculation of spontaneous emission inside such 
a medium would involve the actual computation of the Green function, which for general e(r, uj) is not 
easy. 

A special case of the former theories is the quantum optical description of inhomogeneous systems at 
frequencies where both dispersion and losses can be neglected. Then a description in terms of modes 
is possible, where the mode functions are harmonic solutions of the classical wave equation featuring a 
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position-dependent dielectric "constant" e(r) [£3|. This encompasses the now theoretically and experi- 
mentally very active research area of the so-called photonic crystals p4| , where a periodic modulation 
of the refractive index at the scale of the wavelength of light can drastically modify the mode structure 
compared to vacuum. By increasing the refractive-index contrast, even a photonic bandgap can open 
up, giving rise to a frequency interval for which waves cannot travel in the crystal in any direction, so 
that spontaneous emission would be inhibited completely. Until now, such a bandgap has not been found 
conclusively in the optical regime pEj ], It has been proposed to look for frequencies close to material 
resonances, where refractive indices can be quite substantially higher or lower than 1 pc| . 

Interesting new effects have been predicted for bandgap systems, such as photon-atom bound states 
and non-exponential spontaneous decay at the edges of the gap (27J] . A current debate is whether the 
Weisskopf-Wigner approximation can be used in the calculation of spontaneous emission near an edge 
of a photonic bandgap. This question seems to depend strongly on the analytic or singular behavior of 
the density of states at the edges of the gap, which has recently been calculated for face-centered cubic 
and diamond- like crystal structures [^8| . If near the edge of the bandgap a large part of the modes has 
a cavity-like structure, producing nonzero dwell-times near the emitting atom, then an emitted photon 
has a nonzero probability of being reabsorbed, which would give Rabi-like oscillations of the atomic 
population that are missed in the Weisskopf-Wigner approximation. 

Non-exponential decay can also be caused by the interference of possible decay-channels: for short 
times after the excitation of the atom, a larger frequency interval of the medium states plays a part in 
the decay process than for later times. Ultimately, only the refractive index at the atomic transition 
frequency plays a role, all in concordance with the energy-time uncertainty relation. This interference 
process already happens for spontaneous emission in vacuum. However, when the medium has a strong 
jump in the density of states around the atomic transition frequency, the interference effect will change 
substantially. 

To separate the latter cause of non-exponential decay from the former, it is interesting to consider the 
spontaneous emission inside homogeneous lossy dielectrics with strong and narrow material resonances, 
where the density of states can also change very rapidly. Here all states correspond to simple plane wave 
modes, so that real reabsorption processes do not play a role. In this article, we use the damped-polariton 
model formulated by Huttner and Barnett Jl^,|l^] to study the interference effects of spontaneous emission. 
If absorption is neglected in the damped-polariton model, then we are left with the Hopfield model of 
a dielectric [p|p9|] , which has a frequency bandgap inside which the refractive index is purely imaginary. 
The analogy between this polariton band gap system and photonic crystals was drawn in [p0|| . 

The organization of the paper is as follows: in sectio n [H| we introduce the theory and solve its equations 
of motion using Laplace transformations. In section pUTwe show that the consistency of our solutions 
depends on the validity of a number of velocity sum rules, which are then proved. In section |Tv| , we find 
that for long times all field operators can be expressed in terms of the initial bath operators, and we give 
an interpretation of the result. We also show how to relate the result to phenomenological quantization 



theories. Before we can discuss transient effects of spontaneous emission in section VI, we discuss in 
section |y| the Lorentz oscillator model and the point scattering model. We show how both these models 
can be found from the damped-polariton theory by choosing a suita ble c oupling to the bath. The paper 



ends with a discussion of the results and with conclusions in section VII 



II. THE MODEL AND SOLUTIONS OF THE EQUATIONS OF MOTION 

The damped-polariton theory describes the interaction of light with an absorbing homogeneous 
medium. The coupling of the matter to a frequency continuum is the cause of the light absorption. 
The continuum could be a phonon bath or something else, but for the moment that is not specified: it is 
a collection of harmonic oscillators with a frequency-dependent coupling to the matter fields. Since the 
medium is homogeneous, the dynamics can be separated into a transverse and a longitudinal part. In this 
article we concentrate on the transverse excitations as described by the following Hamiltonian pJLH : 



with 



H — H em + #mat + #bath + Hint, (1) 

H em = J d 3 k hkca j (\,k,t)a(X,k,t), (2) 

tfmat = J d 3 k too fo f (A, k, t)b(X, k, t), (3) 

H hath = J d3k l dwfiw6i(A,k,i)6 w (A,k,t), (4) 
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d J fc 



dej hV(u>) [b(X, k, t) + fet( A , k, t)] [bt (A, k, t) + & W (A, -k, t)] 



J / d 3 fc hu e J^ [o(A,k,t) +a t (A,-k,t)] [&t(A,k,t)-6(A,-k,i)]- 



(5) 



We use the same notations as in [Q. In particular, fee stands for ^Jk 2 c 2 + ui 2 , where the frequency ui c 
equals a/y/pio, with a the coupling constant between field and matter, and p the density. The resonance 
frequency loq of the polarization field is renormalized to ujq, which is the positive- frequency solution of 



LVr, 



uJq + luo / dui V 2 (u})/uj. 



(6) 



The fc-intcgrals in the Hamiltonian arc understood to also denote a summation over the two transverse 
polarization directions labeled by A. The creation and annihilation operators satisfy standard bosonic 
commutation relations. The Heisenberg equations of motion for the bath annihilation operators are: 



MA, k, t) = --V{u) [6(A, k, t) + 6t(A, -k, t)] - iub u {\, k, t), 



(7) 



and similarly for the creation operators. In the following we drop the (A, k)-labels. We solve implicitly 
for the bath variables, as was done in [till in a classical treatment of the model: 



b u {t) = ~V(u) / dt'[b(t') + tf(t')] e -Mt-t') +&w (o) e - 



(8) 



The annihilation operators are defined in terms of the (transverse) physical fields: 

a(t) = 
b(t) = 



2hkc 



kcA(t) - iE{t) 



P 



u> X(t) + -P(t) 
P 



(9) 



and similarly for the creation operators. Here A and E are the vector potential and the electric field, 
X the polarization field and P its canonical conjugate. Insertion of the solution (g) and its Hermitian 
conjugate in the equations of motion gives: 

E(t) = (kcfA{t) + (c c »P(i), 
A(t) = -E(t), 

X(t) = (cj c /a) 2 e P(t) + (u; 2 c /a)eoA(t), 



P(t) = -a 2 ul/{s wl)X(t) + a 2 ti /{2s ul) / At' F(t-t')X(t') - B{t). 

Jo 

In the last equation, the bath operator B{t) is defined as 

B(t) = r dwi V>0 [ b ui (0) 4 + bl (0) 



whereas the the function F in the convolution in Q1Q) is: 



iuj\t 



F(t) 



duii V 2 (uii) sin(wit). 



(10) 



(11) 



(12) 



We get a system of algebraic equations by taking the Laplace transform, which we denote by a bar: 

< 2 Ja \ 



(p 


-l 2 c 2 





1 


P 








-e uj 2 c /a 


P 




a 2 Q 2 /e uj 2 


[1- 





-e Q uj 2 /a 2 



( E(p) \ 
A( P ) 
X{p) 
V P(P) J 



E(0) 
AQO) 
X(0) 
P(0) - B{p) 



(13) 



Through the operator B{p) the bath remains part of the system of equations: this is as far as we can 
"integrate out" the bath variables. 

Now we can determine the dielectric function s(lo), which is a classical quantity, by putting the deter- 
minant of the (4 x 4) coefficient matrix to zero. The determinant gives the dispersion relation 
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D{p) = e(p) P 2 + k 2 c 2 = 0, 



(14) 



with the "Laplace dielectric function" 



+ Q 2 - \w Q F(p) 



(15) 



The function F(p) is the Laplace transform of F(t), which was defined in equation (jl2|)- From this we 
find the dielectric function 



e{ui) — e(— iuj + n) 
with infinitesimal positive rj and 



uj 2 - Co 2 + \CjoF{uj) 



Fiuj) = F(-ioj + n)= dwil/ 2 (wi) 



1 



1 



uj± — uj — in ui\ + uj + in 



(16) 



(17) 



The difference between F{u>) and is denoted by their arguments. The dielectric function satisfies 
the Kramers-Kronig-relations and has the property of a response function that e{— uj*) equals £*iu>). 
It can be shown that it has no poles in the upper half plane, provided that the integral in (0) exists. 
Previous authors |l^,|3l],|32j assumed that the analytical continuation of V 2 (w) to negative frequencies is 
anti-symmetrical in frequency. Then (^) reduces to the dielectric constant in |3l}| , where it was shown 
to be identical to the more complicated expression in H ], 

We combine ( |l3| ) and (|l5|) and write the Laplace fields in terms of the fields at time t = 0, with 
coefficients that are functions of the Laplace dielectric function e(p) and susceptibility x(p) = ^ip) — !• 
For the electric field we find: 



E(p)^D- 1 ip) (pE(Q) + [p 2 xip) 



, 2 2 

K C 



A(0) 



-P 



r p2 



x(p) ^ 1 



XiO) + -p 2 xip)[PiO)-Bip)} 



Q 



(18) 



The other Laplace operators can be found in the same way and are listed in the Appendix. The inverse 
Laplace transform gives the fields at time t in terms of the fields at time t = 0: 



E{t) = M EE it)EiO) + M EA it)AiO) + M EX (t)X{0) + M EP (t)P(0) + B E {t), 
where, for instance, 



M EE it) 



i pioo 

— dpe^ D- 1 ip)p. 



(19) 



(20) 



The operator B E (t) in equation (19) is the contribution of the t — bath operators to the electric field. 



This term will be analyzed in more detail in section IV 



The equal time commutation relations of the field operators are 

[A{\, k, t), -e E(\', -k', t)] = [X(A, k, t),P(X', -k', t)} = ihdw 6 (k - k') 



(21) 



All other inequivalent combinations of operators commute. In particular A and X are independent 
canonical variables. Hence, we have the propert y A, — D] = [A, — EqE], with the displacement field D 
defined as 6qE — aX. With the help of ( |l9| ) and (pi), we can also calculate non-equal time commutators, 
for example: 



[E(X, k, t),E(X', — k', 0)] = M EA (t)[A(X, k, 0), E(X', -k', 0)] 



ih 
so 



Af BA (t)<J A v«S(k-k'). 



(22) 



In principle we have solved the complete time evolution of the field operators. In section III we analyze 
in more detail their short-time behavior, whereas in section [fv| we consider the long-time limit. 
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III. SHORT-TIME LIMIT: SUM RULES 



For fixed fc, the zeroes of the dispersion relation (|lj) are the poles of the integrand in (|2C|). We assume 
that they are simple first-order poles and rewrite the integral (20) as an integral over frequencies u> — ip. 
Then, using contour integration in the lower frequency half plane, we find the coefficients for the electric 
field: 



Mee if) 
M EA {t) 
M EX {t) 
M EP {t) 



= ^Re 

3 



V P,3 V 9,3 -iCljt 
1 



— fcc^^ Im 

3 



UK C sr-^ 

3 



V P-J V 9,3 



Re 



kc 



E 



Im 



"(I 



I - 



n 2 



_P<3 \ c -»Ojt 
n 1 



k 2 c 2 



(23) 



Some details of the calculation and a list of coefficients M mn (t) of other operators can be found in the 
Appendix. In these expressions, the frequencies flj = £lj(k) are the complex-frequency solutions of the 



dispersion relation u> e(w) — fe e = 0. All flj(fc) have a negative imaginary part. Since e(- 



it follows that — f2*(&) is also a solution of the dispersion relation. We can choose to be the 

solution with a positive real part. The summation over j is a summation over all the polariton branches 
of the medium. For each branch, the complex phase velocity is defined as v p j(k) = flj(k)/k and the 
group velocity as v g j(k) — dflj(k)/dk. For convenience, we leave out their explicit fc-dependence in the 
following. 

From equation (|l]|) we can see that the "diagonal" coefficient MEsit) in (p3| ) should have the value 
I at time t = and the "off-diagonal" coefficients M E a(0), M E x(0), etc. should have the value 0. The 
coefficients of the other field operators should also follow this rule. If these constraints are satisfied the 
non-equal time commutators like ( p2| ) get the right equal-time limits as well. The coefficients (^3|) can 
only have the right t = limits, if certain velocity sum rules are satisfied. 

Velocity sum rules can be derived in a systematic way by evaluating the following two types of integrals: 



dw 
do; 



(w + iS) n 
s(u})uj 2 — k 2 c 2 
(to + i5) m X {io) 
e(lo)uj 2 — k 2 c 2 



for n = — 1, 0, 1, 
for m = -1,0, 1,2,3. 



(24) 
(25) 



Here e(co) is an arbitrary dielectric function that satisfies the Kramers-Kronig relations, so it is not 
necessarily of the specific form ([l6|). The integrals can be evaluated using contour integration in the 
complex frequency plane. We can close the contours either in the upper or in the lower half plane. 
Equating the two answers gives a velocity sum rule. In this way one finds for all wavevectors k: 



^Re(u 9j u pJ /c 2 ) = 1, 

3 

E H,,i '''<• •<' ' = L 



(26) 
(27) 



These sum rules can be found from 
obtained before 



|) with n — 1 and n — — 1, respectively. Both relations have been 

The second was coined the Huttner-Barnett sum rule in j33|, because of its 

A second group of sum rules has the 



importance in phcnomcnological quantum theories of dielectrics 
form 



E Im ^ 



v 2g ) 
g,3 u P ,jj 







Vq 



-1,0,1. 



(28) 



The rules with q = —1 and q = 1 follow from ( |25| ) with m — and m = 2, respectively; the case with 
q = follows from (|24|) with n = 0. 

All of these sum rules are independent of any specific form of the dielectric function, as long as it 
satisfies the Kramers-Kronig relations. Other sum rules do depend on the behavior of e(w) for high or 
low frequencies. For example, from fcq) with m = — 1 we find: 
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3 



Re(c\, 3 /vl 3 ) = e(0). (29) 



This sum rule depends on the static limit of the dielectric function. For conductors the dielectric function 
is singular at oj = J35| , but for dielectric functions which can be found from the damped-polariton model, 
e(0) is finite. Two other sum rules can be derived when for high frequencies <x> 2 x(w) approaches a constant 

Aim- 



3 



: Re(%j<4/c 4 ) = 1 + (^lim/fcc) 2 . (30) 

Moreover, if w 2 x(w) + ojf im falls off faster than oj^ 1 , then the integral 

oj 2 \oj 2 y(oj) + oj? 1 

dw 7 r 2 2 ( 31 ) 

, e\oj)oj a — fc z c z 



produces the sum rule ( |28| ) with q = 2. 

Returning now to the time-dependent coefficients ( |23| ) (and the other ones in the Appendix), one finds 
by inspection that one needs all the above sum rules except ( |29| ) to prove that the coefficients have the 
right limits for t = 0. In particular, from equation (|l^) it follows that the fre que ncy 0Jn m as defined above 
exists in the damped-polariton model and equals oj c . Then with ( p7| ) and (|3^) we see that indeed one 
has M E x(0) = in (||. 

It is easy to prove the above sum rules in the following one-resonance model: 

E ( W ) = 1-^L T . (32) 

This e(oj) is real and violates the Kramers-Kronig relations, but it can be considered as a limiting case 
of an acceptable dielectric function. The high-frequency limit of oj 2 x{oj) indeed equals — oj 2 . The two 
sum rules (p6|), ( p7| ) were shown to be valid for this model JlJ] and we want to check (|3(]) as well. The 
dispersion relation is 

oj 4 - {ojI +oj 2 c + k 2 c 2 )oj 2 + k 2 c 2 oj 2 = 0, (33) 

which has two (real) solutions and fi 2 , with sum (oj 2 +oj 2 + k 2 c 2 ) and product k 2 c 2 o> 2 . It follows that 
for all k 



in agreement with (|30|). The other sum rules can also be checked for this simple model. The sum rules 
( psj ) obviously hold, because all group and phase velocities are real in this model. In models that respect 
the Kramers-Kronig-relations, these sum rules are nontrivial. 



IV. LONG-TIME LIMIT 



A. Field and medium operators 



The coefficients MEE{t) etc. in (p 
own characteristic damping time Tj{ 



damp out exponentially in time. Every polariton branch has its 
= l/(Im ttj(k)). After a few times the maximum characteristic 
damping period, with the maximum taken over all branches, the exponentially damped coefficients can 
be neglected. We call this the long-time limit. The s pee d at which it is attained, depends on e(oj) and on 
k. For long times, only the bath operator B^t) in (19) survives, because it has poles on the imaginary 
axis in the complex p-plane: 



B E (t) 



lixiojc 



hojQ 
2^ 



dwi V(oji) 



dp e 



Pi 



p 2 x(p) 



e{p)p 2 



MO) , 

p + ioj\ p — ioj\ 



(35) 



Hence, in the long-time limit, all field operators are functions of the initial bath operators alone. For the 
electric field we find 



£(*)-> .E,(t) 



1 

0J r _ 



2^ 



dwi V(oj\) 



^ 2 xK)M0)e~^ 
e(oJi)oj\ — k 2 c 2 



iujit 



+ 



e*(oji)oj 2 — k 2 c 2 



(36) 
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where the subscript I denotes the long-time limit. The temporal (and spatial) Fourier components of the 
long-time solutions are: 



+ i_ //^o vn^xnuo) 

' {UJ> uo c V 2e s{u)lu 2 - k 2 c 2 ' 



ATM 



uj c V 2e s(lo)uj 2 — k 2 c 2 

1 ltub e V{Lo)(u 2 -fc 2 c 2 )xHM0) 



aoj c V 2 e(cj)o; 2 — fc 2 c 2 



+ ia /fa^" y(a;)(a; 2 - k 2 c 2 - t^VxHMO) . . 

lo% V 2eo £[uj)lo z — fc z cr 

where the superscript + denotes the positive-frequency component of the operator. For future reference 
we also give the long-time limit of the electric field operator as a function of position and time: 



Ej(r,t) 



2(27r)3 £oW 



2 



J dk J2 e A(k) J dwi 



A=l,2 



e(wi)ti; 2 — k 2 c 2 



(38) 



Similar expressions can be given for the other operators. Notice that these long-time solutions indeed 
are solutions of the equations of motion (|io| ) and of the Maxwell equations. The canonical commutation 
relations ( |2l| ) should be preserved in this long-time limit. Also, the non-equal time commutation relations 
like in equation ( p2| ) should be time-translation invariant. The commutation relations can be verified with 
the equality 

^F 2 H |xM| 2 = Im xM = Im = (39) 

which follows from equations (|l6| ) and (0). Since £i(ui) is anti-symmetric in w, all commutators can be 
shown to be proportional to integrals over the whole real frequency axis. Contour integration then leads 
to the required results. 

The solutions found above can be related to those obtained by explicit diagonalization of the full 
Hamiltonian of the model. In p| this diagonalization was carried out by using Fano's technique. In that 
way the field and medium operators were written in terms of the diagonalizing annihilation operators 
(called C(k, u) in H]) and the corresponding creation operators. If one replaces the bath annihilation 
operators b w (k, 0) in the long-time solutions ( |37| ) by the diagonalizing annihilation operators C(k, cj), and 
if one makes similar replacements for the creation operators, the expressions for the field and medium 
operators in are recovered. 

The long-time solutions can be interpreted as follows: when the dielectric medium is prepared in a 
state that is not an eigenstate of the Hamiltonian and if the coupling V(u>) is nonzero for all frequencies, 
then the medium tends to an equilibrium that is determined by the state of the bath. The time it takes 
for this equilibrium to settle down is the time after which the long-time solutions can be used for the 
field operators. So one can always use the long-time solutions in the calculations, unless the medium 
has been specially prepared in a non-equilibrium state a short time before one does the experiment. The 



interpretation of the long-time solution will become clearer in section VI where we calculate spontaneous 
emission. 

In summary, for times long after t = 0, all field operators can be expressed solely in terms of the 
bath operators at time t = 0. The time evolution is governed by the bath Hamiltonian alone. The field 
operators still satisfy Maxwell's equations and the canonical commutation relations. Classical expressions 
for the Maxwell fields would have died exponentially to zero in this long time limit. 



B. Relation with phenomenological theories 



The long-time solutions of the field operators can be related to expressions in phenomenological theories, 
as we will show presently. In phenomenological quantum mechanical theories of homogeneous absorbing 
dielectrics Q - pi] , a noise current density operator J is added to the Maxwell equations in order to 
preserve the field commutation relations: 

VxE + (r,w)=wB+(r,w), (40) 
V x B + (r,w) = -i^ D + (r,^) + ^ J + (r,^)- (41) 
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The displacement field D + in the last equation is defined in terms of the electric field and the dielectric 
function as 



D + (r,w) =e e(w)E+(r,a;). 



(42) 



We write D to stress the difference with the microscopic displacement field D in section [Tl| After taking 
the spatial Fourier transform, and using B + = V x A + and E + = iusA + , so that the first of the Maxwell 
equations is satisfied, one finds from the second equation: 



[ w 2 e{u) - k 2 c 2 1 A + (A,k,w) = J+(X,k,Lu) 



(43) 



The vector potential and all Maxwell fields can be calculated in terms of the noise current density J. 
The canonical commutation relations are preserved, if for the noise current one chooses Mj2(J : 



[ J+(A,k, w ),[J+(A',k>')]t] 



Instead of using the noise current operator, one defines basic bosonic operators 

/(A, k, u) = / n J+(A, k, w), 

so that these operators satisfy simple commutation relations: 

[ /(A, k, «), P (A', k', uj 1 ) ] = 5 xx ,&{k - k')5(cj - u '). 



(44) 



(45) 



(46) 



Now we turn to the long-time solutions of the field operators that we determined in section [V A. The 



long-time solution of the vector potential in ( |37| ) obviously is a solution of the following inhomogeneous 
wave equation: 



[ e(Lup 2 - k 2 c 2 ] A+(A,k,w) = -\/^ n«)wxMMA,k,0). 



(47) 



This kind of equation is well-known in Langevin theories |T^,|l8| : the coupling to a bath gives a damping 
term (here: a complex dielectric constant) in the equations of motion of the system. Besides damping, 
there is an extra term that is neglected classically. This term is the quantum noise operator, which 
features the bath operators at time t = 0. 

The long-time solution ([47]) can justify the phenomenological equation if we identify 



\V(w)x(u)\ 



(48) 



where we used equation (|39|). We see that up to a phase factor, the bath operators fr^A, k, 0) from the 
microscopic theory serve as basic bosonic operators /(A, k, w) in the phenomenological theories. We want 
to stress that the identification (^8|) is only valid in the long-time limit when the medium is in equilibrium 
with the bath. 

In section il| we saw that —£qE is the canonical co njug ate field of A and that [A, —D] gives the canonical 
result as well. Since we can make the identification (pq), the same relations hold in the phenomenological 
theory that was described in this section. But now let us calculate the commutator [A, —D] with D + 
defined as in equation (^) and D~ as its Hermitian conjugate. We can use the long-time solutions, 
because the commutation relations are preserved: 



A(X, k, t), —D(\', -k', t)] = —6 xx ,S(k - k' 



dw 



£ r (u>)Ei(u))uJ 3 

\£{uj)lu 2 - k 2 c 2 \ 2 ' 



(49) 



with e r (ui) the real part of the dielectric constant. The symmetry of the integrand enables us to rewrite 
the right-hand side as an integral over all real frequencies. When using contour integration, one cannot 
replace £*(u>) by e(— u>*), but the analytical continuation to complex frequencies of £*(u>) = e(—oj) must 
be used instead: 



A(X,k,t),-D(X',-k',t)] =iM A A'5(k-k')^Re [e(-ft> PJ% ,/c 2 ] , 



(50) 



where we assumed as before that all poles of the dispersion relation are first-order poles. Note that e(— fij) 
depends on the behavior of the dielectric function in the upper half plane. Contrary to a statement in 
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|2CJ, the commutator does not give the canonical result, because in general there is no sum rule for 
the right-hand side of the equation. In other words, (D — EqE) is canonically independent from E, but 
(D — SqE) is not. The operator (D — D) is proportional to the Langevin noise term in the wave equation 
for the electric field. 

Now let us neglect absorption at all frequencies. Strictly speaking, the limit £i{uj) — * is unphysical 
because it violates the Kramers-Kronig relations, but the limit is sometimes taken for dielectrics that 
show negligible absorption at optical frequencies [{Hps] ]. When s{io) becomes real, the solutions ilj 
become real and in that limit one has e(— flj) — > e(0^) = (c/v p j) 2 . Inserting this in (|50| ) and using the 
Huttner-Barnett sum rule ^ - Re {v g j/v p j) — 1, we immediately find the canonical result for [A, — D\. 
We compare this with the results in |3j|, where the dielectric function is assumed to be real. There 
a phenomenological Lagrangian was introduced and the fields A and — D were correctly identified as a 
canonical pair. The Huttner-Barnett sum rule was invoked to show that their commutator indeed had 
the canonical form. It was concluded that it is misleading that also [A, — sqE] has the canonical form. 
Here we have learnt that this misleading result is not surprising: in the limit of real dielectric constants 
and only then, both [A, — £qE] and [A, —D] can have the canonical form in the same gauge, the reason 
being that D approaches D in that limit. 



V. MODEL DIELECTRIC FUNCTIONS 



Phenomenological theories as discussed in section IV B have expressions for e(io) as input. In practice, 
this input will be the outcome of measurements of the dielectric function. By choosing the appropriate 
microscopic coupling constants and resonance frequencies in the damped-polariton model, one can hope 
to find a given dielectric function, thus providing a connection with phenomenological theories. It was 
argued in pi] ] that the well-known Lorentz oscillator form of the dielectric function could not be found 
from the damped-polariton theory in this way. We shall reconsider this issue below. 

A dielectric function that follows from the damped-polariton Hamiltonian (Q) will have a single reso- 
nance, because there is only one resonance frequency loq in the matter fields. Experimentally, one may 
find more resonances in the e(u>). This should not be used as an objection to the damped-polariton 
model, because in principle one could easily extend the theory with more material resonances. In this 
section, we consider two of these one-resonance models. 



A. The Lorentz oscillator model 



We want to find microscopic coupling constants in the damped-polariton theory so that the resulting 
e(lo) has the following Lorentz oscillator form: 



ELorM = 1 



:,Lor 



■ 2itUK,Q 



(51) 



Here w res is the resonance frequency of the medium and lu c ,Loi- is a frequency that is related to the 
coupling strength between the electromagnetic and the matter field. Identifying e(lo) from equation ( |l6| ) 
with £Lor(^), we find apart from the trivial identification lo c = w C) Lor 



dwiF 2 (cji) 



1 



1 



LOl — LO — IT] LO\ + LO + IT) 



= 4 



LO 



w 



IK + A, 



(52) 



where the frequency shift A is defined such that co^ es = lOq— lo A/2. The coupling V 2 (loi) is fixed by the 
identification of the imaginary parts and for all frequencies it equals V 2 {loi) — Akqloi/ (ttloq). However, if 
we insert this coupling in the equation for the real parts, we find that the frequency shift A is infinitely 
large. Also, the renormalized frequency loq in equation (||) blows up. We can solve this problem by 
introducing a frequency cut-off in the coupling, namely V 2 {lo\) = Ak{loi)lo\ / (ttloq) with 



k{loi) 



kqH 2 



n 2 



(53) 



With this choice one finds u>o = \/lOq + 2ko£1, which clearly has a strong dependence on the cut-off 
frequency. The shift A becomes both finite and frequency-dependent: 



AH 



ttloo 



dcji loik(loi) 



1 



1 



LOl — LO LOl + LO 



ft 

= 4k(u)) — . 

LO() 



(54) 
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The principal value integral can be evaluated by means of contour integration in the complex frequency 
plane. In this way we arrive at the following expression for the dielectric function: 

uj 2 

e(oj) = 1 - 5 — — - — — — : — — . (55) 

0J Z — UJq — 2\1[Kq — K{LJ)\ + 2lUlK(0J) 

We can choose ft arbitrarily high (but finite). Quite unlike u>o, the optical resonance frequency w rcs 
approaches ujq from above, the higher we choose the cut-off, since one has: u> 2 es ~ uo 2 + 2kqlu 2 /VL. Note 
that the dielectric function (p5|) has the right high-frequency limit u 2 x{u) — > — to 2 as required in section 



[II. 



It is well-known that there are two branches of solutions of the dispersion relation when the dielectric 
function is of the form (]5l]): there is an upper and a lower polariton branch. The dielectric function ( |5q ) 
gives rise to another branch: it has a purely imaginary frequency with magnitude of the order of the 
cut-off frequency. This "cut-off branch" has negligible fc-dependence. In fact, the leading fc-dependent 
term for large fl is 2iu> 2 Kok 2 c 2 /f2 4 . Clearly, the group velocity on this branch is practically zero, so that 



the contribution of the cut-off branch to the velocity sum rules of section [II can be neglected. 

We conclude that high cut-off frequencies can be chosen such that in the optical frequency regime 
the dielectric function cannot be discerned from a Lorentz dielectric function with resonance frequency 
u>tes = and damping constant kq- The solutions of the dispersion relation of the upper and the lower 



polariton branch together satisfy the sum rules of section [II 



B. The point scattering model 

In general the dielectric function e(to) describes the propagation of a coherent light beam in a fixed 
direction in an isotropic medium. A complex s{lo) means that there is extinction, which can be caused 
either by scattering or absorption, or both. The dielectric function does not contain information about 
the extinction mechanism. A well-known dielectric medium showing polariton behavior is the dilute gas, 
which can be described as a collection of point dipoles that scatter light independently. If only one type 
of elastic scatterers is present, each having only one resonance, then the dielectric function is given by 
|3l: 



4irc 2 T e n 

Jir e u; 3 / C ' 



where n = N/V is the density of the scatterers (not to be confused with the refractive index n(uj)) 
and T e = e 2 /(4:neom e c 2 ) is the classical electron radius. This dielectric function can also be found if 
one supposes that the medium consists of classical harmonically bound point charges whose motion is 
described by the Abraham-Lorentz equation. The dielectric function (|56|) has the property that the 
corresponding T-matrix t(u>) satisfies the optical theorem, with t(u>) defined as s(u>) = 1 — nt(u>)(c/uj) 2 . 
However, (|56| ) is not a proper response function, since it has a pole near the very large positive imaginary 
frequency 3ic/(2T e ). This can be related to the need for the a-causal phenomenon called pre-acceleration 
to avoid so-called runaway solutions of the Abraham-Lorentz equation |37f| . 

Although we know that in the damped-polariton theory only proper response functions can be 
found, we proceed like in the previous subsection and try to find coupling constants that in the op- 
tical regime give rise to the dielectric function (|56|). Equating with ( fT(| ) we get lu 2 = A-Kc 2 T e n and 
F 2 (wi) = 4r(wi)o;?/(37rd;oc), with 

r W - ^ . (57) 

Here we have inserted a convenient frequency cut-off from the start in order to keep finite the frequency 
u)q and the shift A. Contour integration gives ujq = Uq + \/2T e Q 3 / (3c) , and 

A(w) = —r(tu)n(n 2 + Lu 2 ). (58) 

6UJqC 

The dielectric function has the form 

, x = 1 Airc 2 T e n 

£{UJ> uj 2 -uj 2 + [V2W(n 2 -uj 2 )/{^ 2 c)][T e -T{io)} + liY{io)^/c 1 ' 

In this case, the resonance frequency shifts to frequencies lower than ujq and the shift is larger for 
larger cut-off frequencies. However, since the classical electron radius is so much smaller than an optical 
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wavelength, it is very well possible to choose a cut-off frequency such that ojq <C £1 <C c/T e . Then 
for optical frequencies, the dielectric function (p9|) is of the form (p6|). Note that for high frequencies 



w x( w ) ~~ * ~ w c f° r the dielectric function (p9[), but not for (p6[). 

Again, the frequency cut-off introduces a cut-off branch. In Fig. |l| we plot the real parts of the three 
solutions fij(fc) of the dispersion relation. As a measure of the damping, we introduce n which is given 
by T c luq I (3c) . For purpose of presentation, the numerical values of both lu c and k were chosen artificially 
large for a dilute gas. The frequencies on the cut-off branch are of the same magnitude as the cut-off 
frequency O, much higher than the optical regime. The imaginary parts of the upper and lower polariton 
branches are plotted in Fig. ^|. The imaginary part of the cut-off branch is large negative and practically 
constant for parameters as given in Fig. [l|. Again, since the group velocity on the cut-off branch is 
practically zero, the upper and lower polariton branches together satisfy the sum rules of section III. In 
particular, Fig. || illustrates that the upper and lower polariton group velocities v 9yU and v g> i satisfy the 
sum rule Im {v 3lU + v g> i) = 0. 




kc/co 

FIG. 1. Real parts of the three solutions flj of the dispersion relation with e(w) as in (fs^). Numerical values 
of the parameters: f2 = ltkjo, w c = 0.5uio and k — O.Ol^o- The solid line is the lower polariton branch, the upper 
polariton branch is dashed and the cut-off branch is dotted. 



0.000 




-0.010 



FIG. 2. Imaginary parts of the lower (solid line) and upper (dashed line) polariton solutions Qj of the dispersion 
relation with e(w) as in (|5^). Numerical values of the parameters as in the previous plot. Not shown is the 
imaginary part of the cut-off branch, which is also negative and about a thousand times larger in magnitude. 



The cut-off, which was necessary to produce the dielectric function in the damped-polariton theory, 
neatly removes the pre-acceleration behavior associated with a pole in the upper halfplane and leads to 
a good response function. The form of the coupling V(u>) given above (|5^) has the following physical 
interpretation. By equating the damped-polariton dielectric function with (|56|), we assumed that the 
dilute gas can be described as a homogeneous dielectric. The light scattering by the gas molecules can be 
accounted for by an absorptive coupling to the free electromagnetic field, as long as only single scattering 
of light is relevant. Then scattered light is lost for propagation in the original direction. If the matter-bath 
coupling is dipole coupling, then for optical frequencies the product V 2 {loi)/uji should be proportional 
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to the density of states of the electromagnetic field, which goes quadratically in frequency. This is indeed 
the case. 



VI. SPONTANEOUS EMISSION 



The spontaneous emission rate in principle is a time-dependent quantity. In this section we investigate 
the transient dynamics of the spontaneous emission rate of a guest atom in an absorbing medium, when 
the transition frequency of the guest atom is close to a material resonance of the medium. We show 
how our results relate to previous treatments of spontaneous emission in absorbing dielectrics, where 
Fermi's Golden Rule was used to show that the time-independent (equilibrium) value for the spontaneous 
emission rate equals roRe[n(wA)] [[ll]jl6). Recently, local field effects have been included in quantum 
electrodynamical formulations of the problem P,p|,p| 1 1 36|] , but we shall not focus on them in this paper. 

We model the guest atom as a two- level atom with ground state \g) and excited state |e) and Hamilto- 
nian Ha = fwA\e) (e|. The medium (with fields and bath included) is described by the damped-polariton 
model, with Hamiltonian Hm given by ([!]). The total Hamiltonian is H = Hq + V, with Ho — Hm + Ha 
and V — —fiA ' E(ryi), the dipole interaction between the atom and the medium; \xa is the atomic dipole 
moment operator and E(r^) is the electric field operator at the position va of the atom. 

Suppose that the damped-polariton system is prepared at time in a state described by a density 
matrix pm(0). We do not assume that pm(0) commutes with Hm, nor that it factorizes into a product 
of a density operator for the bath and a density operator for the undamped-polariton system (as is often 
assumed for convenience p9[). At time to > we bring the guest atom in its excited state, and couple it 
to the damped-polariton system. Using perturbation theory, one can calculate jl8| the time-dependent 
probability that the guest atom has emitted a photon at time t > to- We define the derivative of this 
quantity as the instantaneous spontaneous emission rate T(t). It is given as 

T(t) = ^Re^ dt' e^-^TTlpMiO)^ ■'Eir A ,t)ti--E(T A ,t')], (60) 

where fj, is now the dipole transition matrix element of the guest atom. 

If the guest atom is excited a long time after the initial preparation of the medium, all transient effects 
in the electric field have damped out. Hence, the field may be replaced by its long-time limit Ei(rA,t), 
which is given in (38). Since E; depends only on the bath operators at t — 0, we may write ( |60|) in the 
form: 



r(t) = ^Re jf dt' e <WA <*-*''Th«thbr e d(0)M-Ei(rA,*)/i-Ei(rA,* / )]. 



(61) 



Here p re d is the reduced density matrix obtained by tracing out the electromagnetic and material degrees 
of freedom: p rc d(0) = Tr em ma t/9M(0). For the special case that the initial density matrix pM (0) factorizes, 
the reduced density matrix is the bath density matrix pbath(O) at t = 0. In general, the initial state of 
the electromagnetic and material degrees of freedom at t = does not play a role in the emission rate. 

Spontaneous emission in its pure form arises if the reduced density matrix describes the ground state 
of the bath. Let us assume this is indeed the case. Upon inserting ([38]) in (|6l]) we can perform the 
t'-integral, the integrals over the wavevector and the summations over the polarization directions. This 
leads to 

r W = T~2i 3 Re / duv n(uj) , 62 

3ir z heoc- i J Q uo — uoa 



with n(u>) = \J e(w) the complex refractive index. 

For times (t — to) that are large enough, one may replace sin[(w — u>A)(t — to)]/ (uj — uja) by ir5(u> — u>a)- 
However, the time scale at which this replacement is valid, depends on the resonance structure of the 
refractive index n(u>). Since we want to study just this time scale, we will not make the replacement. 
To evaluate the integral we multiply the integrand by a convergence factor ft 4 / (£! 4 + us 4 ), with fl 3> oja- 
The specific choice of the cut-off frequency £1 will only affect T(t) at time differences t — to much smaller 
than a single optical cycle. We need to use a high-frequency cut-off at this point, because the dipole 
approximation is incorrect for high frequencies. 

For the dielectric function, we take the Lorentz oscillator form (p5|), and we choose the cut-off frequency 
in that model to be identical to the one inserted in (^). In Fig. |pwe give the real part of the refractive 
index which clearly changes rapidly near to = ujo- It is a familiar figure and it shows that the refractive 
index does not change much while increasing the cut-off frequency Q from IOuq to infinity. The density 
of radiative modes around the material resonance is proportional to uj 2 Ke[n(uj)]. 
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co/co 

FIG. 3. Real part of the refractive index in the Lorentz oscillator model, when e(w) is given by equation ( |55| ) 
with parameters ko = O.Olwo, = 0.5o;o, Q — Wu>o. 

With this model for the dielectric function and the parameters as in Fig. ^, we calculated T(t) in the case 
that the transition frequency u>a exactly equals ljq. Since the integrand in ([32]) is rapidly fluctuating, it is 
expedient to use complex countour deformation to evaluate the integral. We add an infinitesimal positive 
imaginary part to the denominator and split the sine into two complex exponentials. The contour of the 
integral with exp[i(w — u>A.){t — to)] in the integrand is deformed towards the positive imaginary axis. The 
contribution from the pole arising from the convergence factor can be neglected at time scales t— to ^S> o;7 ■ 
Likewise, the integration contour of the integral with exp[— i(co — LOA)(t — to)] is deformed towards the 
negative imaginary axis. Again, the pole contribution from the convergence factor is negligible. Further 
contributions, which cannot be neglected, arise from the branch cuts of n(uj) = \/e(iJ) and from the pole 
at loa- The latter contribution is easily evaluated and yields the equilibrium value r(oo) = To Re n(u)A)- 
In contrast, the branch cuts yield time-dependent contributions to T(t). For large f2 they are situated at 
uji = —ino + y^uiQ ~ K o an d ^2 = — ii^o + ^/o^q + cj^ — Kq. Around u>\ and 0J2, we can approximate the 
dielectric function by 



e{uj 2 



-2iSe 



The branch cut at lo\ gives the following contribution to the spontaneous emission rate: 



-Re 



-i7r/4-re (*-to)-«(\/ w o _K o _CJ ^)(* _t o) 



Jit) 



(63) 
(64) 

(65) 



where J(t) is defined as: 

e-^-^i^Jf 1 ^M(A + k )] 3 



J(t) 



dX _ 



lu a - i(X + /to)] 



2u 2 c y/ul - n\ + zA(A 2 + Amp 1 + J\ - 4 k 2 ) 
4(a; 2 - K 2 ) + A 2 



l 1/2 



(66) 



The branch cut around loi gives a similar contribution. 

The integrals arising from the branch cuts and from the imaginary axis can easily be evaluated numer- 
ically, since their integrands are no longer rapidly fluctuating. The result is the solid line in Fig. |[ We 
see that the spontaneous emission rate builds up until it finally reaches the time-independent equilibrium 
value Tq Re n[u>A)- 
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FIG. 4. Normalized spontaneous emission rate F(t)/[To Re tl{uja)] in an absorbing dielectric as a function of 
time, when the transition frequency lua is equal to c^o- Choice of parameters in the Lorentz oscillator model: 
no = O.Olwo, cOc = 0.5ojo, ^ = oo. The solid line is the exact result for ( p^ ) and the dashed line is the approximate 
expression (p7[). 

The dashed line in Fig. |] is an analytical approximation for T(t), which captures the main features of 
the time dependence, at least qualitatively. It is derived by retaining only the contribution (^5|) in the 
time dependent part of r(t), as this is dominant for large t. Moreover, we approximate J(i) by the first 
term in its asymptotic expansion for large t — to. In this way we arrive at the following approximate 
expression for T(t): 



T{t) ~ r Re 



n(Lu A ) 



uj c {ujI 



-in /4-K (t-*o)-i[V '"o _K o _w a](*-*o) 



2iruj\{yjujl - nl - UJ A - IK ) 



(67) 



As explained, this approximation contains only the contribution from the branch cut at w\\ the branch 
cut at W2 gives a faster decaying term, which goes like e~ K °(*~*°) /(t — to) 3 / 2 . The contributions from the 
integrals along the imaginary axis decay even faster. 

It can be seen from equation ([57]) that the amplitude of the time-dependent part of T(t) falls off as 



-K (t- 



*°'/(t — to) 1 / 2 and also, that the amplitude of the extra term is largest around resonance, when 
. Away from resonance oscillations with frequency \/ uJq ~ K o ~ U A are present. Fig. ^ 



o 



shows the on-resonance case when the time-dependent term shows no oscillations, but has relatively large 
amplitude. 

The main result of the present discussion is the time-dependence of the spontaneous emission rate. 
The time-independent value is not reached instantaneously, but at a time scale that is governed by the 
resonance characteristics of the medium. In fact, the smaller the resonance width Ko, the longer it 
takes to reach the time-independent value. Typically, it takes u>o/ K o optical cycles, as follows from the 
exponential e - K o(t-t ) m approximate expression (|6"ij). For narrow resonances with loq/ko large, the 
transient dynamics may take a substantial amount of time. 



VII. DISCUSSION AND CONCLUSIONS 



We have solved the equations of motion for the field operators in the damped polariton model using 
Laplace transformations. The solutions of the field and medium operators are the sum of a transient 
and a permanent part. The latter are expressed solely in terms of the initial bath operators. Long after 
the initial time all field and medium operators are functions of the bath operators alone, provided the 
coupling to the bath is nonzero for all frequencies. The long-time solutions satisfy quantum Langevin 
equations in which the initial bath operators figure as the quantum noise source. The same continuum 
that produces the absorption also forms the noise source that keeps the commutation relations in order. 
This is conceptually simpler than expressing the quantum Langevin noise in terms of the creation and 
annihilation operators that diagonalize the total Hamiltonian of the damped-polariton model fl4|| . 

The effects of the initial state of the field and medium variables on the expectation values at a later 
time are noticeable only during a short period that is determined by the characteristic relaxation times 
of the damped polariton modes. Once these transient effects have died out the expectation values are 
determined by the reduced density matrix which follows from the full density matrix at the initial time 
by taking the trace over the degrees of freedom of field and matter (without bath). If the full density 
matrix at the initial time factorizes, the reduced density matrix equals the initial bath density matrix. 
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The method of long-time solutions can be used for other dissipative quantum systems as well. For 
models in which the Hamiltonian can be diagonalized completely, it is an alternative to the Fano diag- 
onalization technique. The latter can be quite complicated |l4| , f40| ], whereas our long-time solutions are 
found after the simple inversion of a 4 x 4 matrix, as one sees from section [n] and |y|. More generally, the 
long-time method may be useful for dissipative systems with a bilinear coupling to a harmonic oscillator 
bath whose dynamics can be integrated out. 

We employed the method of long-time solutions to study transient effects in a medium described by 
a Lorentz oscillator dielectric function. This dielectric function (and that of the point-scattering model 
as well) can be derived from the damped-polariton model by taking a suitable bath coupling. Although 
a cut-off procedure turns out to be indispensable, the essential physics in the optical regime can be 
represented adequately in this way. Once the connection with the damped-polariton model has been 
established, spontaneous emission processes by a guest atom in a Lorentz oscillator dielectric can be 
investigated by means of the long-time method. Although transient effects due to the initial preparation 
of the dielectric have damped out after a few medium relaxation periods, transient behaviour of a different 
type shows up in the initial stages of the decay process. This transient behaviour, which is related to the 
preparation of the guest atom in its excited state, leads to a non-exponential decay - or in other words 
to a time-dependent spontaneous emission rate - if the atomic transition frequency is near a resonance of 
the dielectric. The non-exponential dynamics takes place at time scales that are inversely proportional 
to the width of the resonance. As we have shown, the characteristics of the time-dependent decay rate 
can be captured in an analytic asymptotic expression of which the qualitative features are corroborated 
by numerical methods. 
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APPENDIX: LAPLACE OPERATORS AND TIME-DEPENDENT COEFFICIENTS 



In section || the electric field operator E(p) was given in terms of the operators at t = 0. Here we 
give the analogous expressions for the other Laplace operators. Furthermore, we show how to evaluate 
the time-dependent coefficients M mn (t) like in (|2^) for the electric field operator. Finally, we list the 
expressions for the coefficients of the other operators. 

The expression for E(p) in (111) has the following analogous expressions for the other Laplace operators: 



A(p) 



Pip) 



Dip) 
a 



X(j>) = 



-E{0)+ P A(0) 



rxO) - 1 



X(0)--px(p)[P(0)-B(p)] , 

a 



D(p) la a 

^ + 4 + A PXip)X i0) + ^i P 2 + k 2 c 2 )x(p){P(0) - B(p)] ) . 



(Al) 



(A2) 



D{p) 



+P 







A- 


-a 


a 2 


y 


£o 




fk 2 c 2 


u + 



1 



-,xip) - 1 



E(0) + ap 



2^,2 



k C 



1 ) x(p)[P(0) ~ B(p)} 



^Xip) 1 



1)X(0) 



A(0) 



(A3) 



If we now apply the inverse Laplace transformation to these expressions, we find the full time dependence 
of the operators A, X and P. The inverse Laplace transformation is a contour integration over the 
Bromwich contour that includes the whole imaginary p-axis. After transforming to frequency variables, 
the contour includes poles from D^ 1 {p), which are in the lower halfplane, and moreover poles on the real 
frequency axis arising from B(p). The latter are important in the calculation of the long-time solutions of 
the operators in section IV. However, in the calculation of the coefficients M mn (t), which we will discuss 
here, they play no role. Let us consider as an example the coefficient Mae^): 
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-200+77 
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2ni 
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4nkc 



duj- 



(IhJ 



n(uj)uj — kc n(u>)u) + kc 



e(u))u) 2 — k 2 c 2 

iuit \ ^ 

kc 



(A4) 



Note that M AE (t) is exponentially damped because all flj in the exponentials have negative imaginary 
parts. The other coefficients can be calculated in a similar way. The results are: 



M AA (t) = M EE {t), 

, x oik 
M AX {t) 



Mxx(t) 



M X p{t) 



M PX {t) 



uj 2 e c 4 



- in ^ 9li (l 



p-j 

^2 



/c 2 c 2 ' 



Map(«) = -Y) Re 
a * — ' 

j 

M XE {t)=E M A p(t), 
M XA (t) =e M E p(t), 



-iQ,jt v g,i ( v p,j _ c 



Soke 



5> 



-iQ.it u g,j 



c 



J P,3 
C 



W P,3 
C 



P,J 
r 2 



M PE (t) =e M AX (t), 
M PA (t) =e M EX (t), 



J2 Im 



P,3 
n 2 



k 2 c 2 



M P p(t) = M xx (t). 



k 2 c 2 



(A5) 
(A6) 

(A7) 

(A8) 
(A9) 

(A10) 

(All) 

(A12) 
(A13) 

(A14) 

(A15) 



With the sum rules discussed in section III, one can see that the "diagonal" coefficients in this list equal 
1 at time t = 0, whereas the other coefficients have the initial value 0. 
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